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ABSTRACT 

We present a detailed analysis of the upper main sequence of the ~1.3-Gyr 
old open cluster Trump ler 20. High accuracy BVI photometry combined with 
the Very Large Telescope/FLAMES medium-resolution spectroscopy of 954 stars 
is essential to understanding the unusual appearance of the color-magnitude di- 
agram (CMD), initially suggesting multiple populations in Trumpler 20. We 
show that differential reddening is a dominant contributor to the apparent split- 
ting/widening of the main-sequence turnoff region. At its extreme, the excess 
differential reddening reaches /\{B — V) ~ 0.1 while the adopted minimum red- 
dening for the cluster is E{B — V) = 0.36. A unique sample of measured projected 
rotational velocities indicates that stellar rotation is high near the main-sequence 
turnoff, reaching vsini ~180 km s~^. By dividing the upper main-sequence stars 
into equal groups of slow and fast rotators, we find that fast rotators have a 
marginal blueshift of 6{V — /) ~ —0.01, corresponding to a difference in the me- 
dian vsini of ~60 km s~^ between these subsamples. We conclude that stellar 
rotation has an insignificant effect on the morphology of the upper main sequence 
of this intermediate-age open cluster. Trumpler 20 appears to contain a single 
coeval population of stars but there is evidence that the red clump is extended. 

Subject headings: open clusters and associations: general — open clusters and 
associations: individual: (Trumpler 20) — stars: rotation — techniques: spec- 
troscopic 



iThis is woes paper 49 of the WIYN Open Cluster Study. 
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Introduction 



A ~1.3-Gyr old Galactic open cl uster, Trumpl e r 20 ( hereafter Tr 20; d = 301?5 b = +2?2 
was first correctly characterized by iPlatais et al.l (120081) . Its bas i c pro pertie s such as age, 



distan ce, and reddening were later confirmed by ISeleznev et al.l (120101 ) and ICarraro et al. 



(l2010l ). These studies also indicate that the red clump of Tr 20 has a peculiar morphology; 
possibly, it is dual. In addition, the upper main sequence appears to show an enlarged 
color spread that is not normally seen in the color-magnitude diagram (CMD) for open 
clusters in this age group. This observational evidence hints at the now-popular conjecture 
of mu ltiple stellar popula t ions, suggested for a number of the Magella nic Cloud (MC) clusters 



(e.g., iMackev et al.ll2008l : Milone et al.ll2009l : iGoudfrooii et al.ll2009f l. Recent studies of rich 
intermediate-age (1-2 Gyr) star clusters in the Large Magellanic Cloud (LMC) indicate that 
their extended main-sequence turnoff ( MSTO) regions can be described by continuo us star 
formation, lasting ~0.3 Gyr or longer (JGoudfrooij et al.l l201ll : iRubele et al.ll201l[ ). This 
long period of star formation, however, seems to be at odds with the fact that none of 
the younger clusters are known to h ave such a trait. Alterna tive interpretations include, 
for exam ple, stellar rotati on effects ( Bastian fc de Mink 12009) and int eractions in binary 
systems (JYang et al.ll201ll ). However, according to iGirardi et al.l ( 1201 if ), stellar rotation is 
not a likely physical process to produce extended MSTO. Although extended star formation 
is generally accepted as the main factor in generating atypical CMDs for a number of MC 
star clusters, this scenario is based solely on interpretation of the photometric data. An 
exception is the analys is of chemical abundances in four intermediate-age LMC clusters 



(Mucciarelli et al. 



20081 



h. Among them, NGC 1783 is a star cluster with an extended MSTO 



(JGoudfrooij et al.ll201ll ): however, it lacks convincing anticorrelation in the Na-0 and Mg- 
Al abundances, obs erved in the Galactic globular clusters with multiple populations (e.g. 



Gratton et al.ll201ll ). This might be an indication that younger star clusters with extended 
MSTO may not experience self-enrichment. 

The open cluster Tr 20 is important for two reasons: (1) if proven, it would be the 
first known open cluster with multiple stellar populations in the Milky Way; (2) the relative 
proximity of Tr 20 {d = 3.3 kpc) permits high-resolution spectroscopy of its upper main 
sequence. This provides the crucial observational data on projected rotational velocities, 
vsini, which are scarce for the main-sequence stars of 1-2 Gyr old star clusters. Such 
data can be used to test theoretical predictions on the significance of stellar rotation on 
morphology of CMD. In this Letter, we discuss the results of high- accuracy photometry of 
Tr 20 combined with measured radial velocities and v sin i for the entire upper part of the 
cluster's CMD. Detailed results of this study will be published elsewhere (I. Platais et al., 
in preparation). 
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2. Observational data 

2.1. BVI Photometry 

The photometric BVI CCD observations of Tr 20 were obtained on three observing runs 
in 2008-2009 at the CTIO 1.0 m telescope, operated by the Small and Moderate Aperture 
Research Telescope System (SMARTS) consortium. The instrumental setup included an 
STA 4064 X 4064 CCD featuring 15 /im pixels and covering a field of view of 20' x 20'. 
Exposure times for Tr 20 ranged from 10 to 600 s. Each night a large number of Landolt's 
standards (57-112) were frequently observed in order to assure the best possible calibration 
of instrumental magnitudes based upon empirical and field-variable point-spread function. 
Fringing in the sky background was substantial in the I filter and was removed accordingly. 
The observed spread of extinction, color, and zeropoint terms in the transformations to 
the standard system indicates excellent photometric quality of the sky. The combined BVI 
photometry for 96,156 stars covers a 43' x 43' area down to \^ ~ 22 with formal accuracies 
reaching 0.003 mag for optimally exposed stars. 



2.2. VLT/FLAMES/GIRAFFE Spectroscopy 

Spectroscopic observational of Tr 20 were conducted in 2009 April-May using the Very 
Large Telescope (VLT) /FLAMES/GIRAFFE multi-object spectrograph in its MEDUSA 
mode, which can allocate up to 132 fibers over a circular field-of-view with a diameter of 
25'. The chosen high-resolution mode H525.8B offers resolution R ~ 25, 900 at the central 
wavelength of A = 5258 A. Our sample of targets contains phot ometric clust e r mem bers 



selected from the cluster's CMD over a field of 20' x 20' from iPlatais et al.l ((20081). It 
includes the entire red giant branch (RGB), all possible subgiants, and the complete upper 
main sequence down to V^ = 17.0, equivalent to My ~ +3.3 (F2-F5V spectral type). The 
CMD of the observed stars is discussed in Section 3. A total of 954 stars were observed using 
nine fiber configurations, each taken with 1 hr of exposure time that resulted in signal-to- 
noise ratio (S/N) of 15-70 for the entire sample. Each target star was observed only once. 
Our goal was to maximize the number of main-sequence cluster members, while ignoring 
potential binaries and concentrating on likely single stars with radial velocities consistent 
with cluster membership. 



The data were reduced using the GIRAFFE BaseLine Data Reduction Software (JBlecha et al 



^Based on observations collected at Paranal Observatory, ESO, program ID 083.D-0671(A). 
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2000l ). which follows through the usual steps of spectroscopic reductions. In addition, we 
modeled the sky and subtracted it from the spectra of main-sequence stars. 



2.2.1. Cross-correlating with Empirical Templates 

Radial velocities V^ and projected rotational velocities v sini were first derived by cross- 
correlating the spectrum of each star with box- shaped- line digital templates, constructed 
from the high- re solution and high-S/ N spectra (FO V and KOIII, accordingly) of slowly ro- 
tating stars (e.g., 



Baranne et al. 



19961 ). The cross-correlation functi on (CCF) was th en fitted 



with a Gaussian whose parameters served in deriving K and vsini (JMelo et al.ll200ll ). These 
techniques work well for RGB stars, which have intrinsically low vsmi values, and for the 
main-sequence stars with vsuii < 80 km s~^. There is a hint that with increasing usinz the 
radial velocities become biased, and for vsmi > 100 km s~^ this method often fails because 
the narrow noise peaks, sig nificantly enhanced by relatively low S/N of our spectra, start 



to dominate the CCF (e.g., iGalland et al.ll2005l ). These shortcomings of the fixed-template 



techniques and the desire to measure higher v sin i cases prompted us to undertake a different 
approach. 



2.2.2. Cross-correlating with Synthetic Spectra 

To determine stellar atmospheric parameters and radial velocities, we used the Harvard- 
Smithsonian Center for Astrophysics (CfA) library of synthetic spectra, which covers the 
entire wavelength range of the GIRAFFE spectra. The present version of th e library was 
calculated by John Laird for a grid of Kurucz model atmospheres (lKurucal2005l ). using a line 
list developed by Jon Morse. We broadened the library to match the actual instrumental 
resolution of the observed spectra, which we determined t o be FWHM=13.3 km s~^ by 
fitting the width of ThAr lines. Similar to the analysis in ICarter et al.l ( l201ll ). we cross- 
correlated each star against the grid of synthetic spectra with the IRAF task XCSAO and, 
using the normalized CCF peak heights at each grid point, we interpolated to the stellar 
parameters yielding the maximum correlation value. We assign errors of half the grid spacing 
— 125 K in Teff, 0.5 - 10 km s~^ in v sini (with smaller errors for slowly rotating stars), and 
0.25 dex in \ogg and [m/H]. Spectroscopic determinations of Teg, \ogg, and [m/H] can be 
correlated, potentially causing systematic errors larger than the nominal errors quoted here, 
but constraining [m/H] to equal the cluster metallicity can help to break that degeneracy. 
We used the red clump cluster members (see Section 3) to determine [m/H]= —0.08 and 
then performed the interpolation again with the metallicity fixed. We also note, as hinted 
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with the other techniques, that the radial velocities are redshifted in the amount of about 
+0.02 km s~^ per unit v sin i value. 



3. Cluster membership and color— magnitude diagram 

The C MP of Tr 20 shows a h eavy contamir i ation by field stars, mostl y from the Carina 
spiral arm (jSeleznev et al.l l2010l : ICarraro et al.l l2010t iPlatais et al.ll2008[ ). Our sample of 
nearly a thousand photometrically selected probable cluster members can now be cleaned 
up by using radial velocities. First, we examined radial velocities (see Section 2.2.1) of the 
red clump stars to derive the mean cluster velocity. A sample of 68 red clump stars yields 
(Vr) = —40.60 ± 0.12 km s~^ and is our adopted mean radial velocity for Tr 20. A similar 
estimate from radial velocities based upon synthetic templates gives (K) = —40.67 ± 0.12 
km s~^. Second, we used vsini and Vr to select the probable kinematic cluster members 
(Figure 1). Our final sample of kinematic members of Tr 20 consists of 471 stars with a 
caveat that among them ~50-100 stars (10%-20%), estimated by two different statistical 
methods (ad hoc analysis of v sin i distribution and a formal cluster membership calculation 
in the bins of vsini), can still be field stars, predominantly at vsini < 70 km s~^. 

The CMD of our kinematic cluster-member sample (Figure 2) shows some intriguing 
features. In order to represent and interpret them, here we chose a V,V — I CMD. First, 
the MSTO region appears to be split into two parts: a narrow but sparse blue sequence and 
a bulky redder one. We note that a split of the upper main sequence in Tr 20 is unusual 



compared to the CMDs of selected LMC star clusters (e.g.. lGoudfrooij et al.ll201ll ). which all 
show a smooth and progressively widening MSTO toward the higher luminosities. Second, 
in the sense of a sequence, the RGB stars are difficult to identify. Essentially, there are only 
two distinct "blobs" of red stars. The most conspicuous one at V^ ~ 14.7 represents the red 
clump, which appears to be significantly elongated along both axes. The second, at l^ ~ 13.8, 
represents the RGB clump, which is well matched by the isochrones (see Figure 4). Prior 
to interpreting these features, we first probed the cluster for signs of differential reddening 
(DR). 



3.1. Differential Reddening 

In CMDs of star clusters older than ~1 Gyr, DR aff ects mostly the mo rphology of the 
MSTO region and the RGB including the red clump (e.g.. |Platais et al.ll201ll ). Our approach 
in estimating DR is based on the assumption that there is an empirical blue envelope of the 
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main sequence, indicating cluster stars with the lowest reddening as shown by the curve 
in Figure 2. This envelope conforms with the adopted youngest isochrone to the cluster's 
CMD. Next, we partially de-reddened each main-sequence star (for Tr 20, down to \^ ~ 17) 
so that they lie exactly on the blue envelope. A zero DR is adopted for a few probable cluster 
members bluer than this envelope. Thus, in our notation, DR is never negative. To mitigate 
the potential correlation between DR and the effects of stellar rotation, we used only the 
stars with vsini < 95 km s~^ (a total of 237). A raw map of estimated DR (Figure 3) 
shows a distinct window of lower reddenin g in th e NW corner of Tr 20, consistent with the 
Galactic reddening maps of ISchlegel et al.l (Il998[ ). The amount of DR is then smoothed on 
a grid of tangential coordinates at 1' steps (Figure 3). At each grid point, the median DR 
is calculated from the nearest few measurements. The average size of a constant reddening 
patch {d ~ 2f5, equivalent to the eight nearest measures) was estimated by minimizing the 
color dispersion in the DR-corrected upper main sequence. After correcting for DR, there is 
a substantial improvement in the tightness of all the features in the CMD (Figure 2). DR 
is clearly the main reason for the unusual appearance of the observed CMD for Tr 20. The 
following isochrone fits to our photometric data indicate that the lowest reddening toward 
Tr 20 is E{B — V) = 0.36 for stars with {B — V)q = +0.3, while the maximum amount of 
DR reache s ^0.1 mag. This range of E(B — V) i s consistent with the earlier estimates of 
reddening (jPlatais et al.ll2008l : ICarraro et al.ll2010l ). We note that the other effects such as 
unresolved or interacting binaries and rotation of stars may partially mimic the appearance 
of DR. 



3.2. Isochrone Fitting 



We used a set of Padova isochrones (JMarigo et al.ll2008l ) generated for the metal abun- 
dance oi Z = 0.015. The i sochroii e s wer e red dened using t h e pre scription of deriving a 
variable E{B — V) given in iFernid (119631 ) and iPlatais et al.l (120081 ). In order to achieve 
a reasonable match of isoc hrones with the features of RGB, the parameter r] defined as 
rj=EB-v{Sp T)/Eb-v(BO) jFernielll963h should be revised to ?] = 0.945 0.17(B - V)n. 



Reddening in BV bandpasses was translated to VI using the coefficients from lSchlegel et al. 



(119981 ). Our best fit is achieved with a 1.3 Gyr isochrone and the true distance modulus 
Vo — M„=12.60 mag. The morphology of the DR-corrected red clump is unusual — it is 
substantially extended along the magnitude axis {AV ~ 0.5), having its highest density 
at V^ = 14.4. There is a certain similarity to a double red clump of the SMC star cluster 
NGC 419; however, the lack of multiple populations in Tr 20 and, hence, the absence of the 
required ^0.1 Gyr s pan of ages appears to weaken the hypothesis of non-degenerate He cores 



(IGirardi et all 120091 ) 



We also used the isochrones from the Dartmouth Stellar Evolution Database ( iDotter et al 



20081 ). A noted 0.2 Gyr age difference from the fit with Padova isochrones (Figure 4) is likely 
due to differences in the treatment of convective core overshoot and updated nuclear reaction 
rates, especially for the CNO cycle. 



Stellar rotation 



Roxburgh fc Strittmatterl (119651 ) were first to predict the effect of stellar rotation on 
the color and magnitud e of a star, claiming that s tellar rotation is essentially tantamount 
to reddening. Recently, iBastian &: de MinkI (|2009[ ) suggested that stellar rotation in stars 
with masses in the range 1.2-1.7M0 can mimic the effect of multiple populations in star 
clu sters by shifti n g the se stars to redder colors in the CMD. This proposition was contested 
by iGirardi et al.l ( l201ll ) on the grounds of evolutionary tracks for non-rotating and rotating 
stars. These authors conclude that stellar models with rotation appear to produce a modest 
blueshift in the CMD. Clearly, our data can provide a crucial test of these predictions. 

First, from the DR-corrected BV CMD we selected a sample of 168 upper main-sequence 
stars down to ^ ~ 16.1 — the part of a CMD containing a relatively high fraction of fast 
rotators. Second, the VI CMD of this new sample was rotated around the point V = 16.1, 
V — I = 0.75 in such a way that the magnitude axis becomes approximately parallel to the 
selected short segment of a 1.3 Gyr isochrone. The rms of the color distribution for the 
modified DR-corrected l^ — / of these stars is 0.04 mag. The modified colors served as an 
argument in constructing the histogr ams of fast and slow rota tors, divided into equal parts 
at vsini = 90.5 km s~^. Similar to iGoudfrooij et al.l (120091 ) . we used the nonparametric 
Epanechnikov-kernel probability density function to smooth the histograms, normalized rel- 
ative to the maximum counts of slow rotators. The resulting color distributions are shown 
for two cases: using VI photometry "as is" and correcting it for DR (Figure 5). Our data 
show that at 1.3 Gyr and in the mass range of 1.5-1.8M0 stellar rotation produces a slight 
to zero blueshift of stars in the DR-corrected CMD: 5{B -V) = 0.0, 5{V - I) = -0.01, and 
6{B — I) = —0.01. All median color differences are between the fast and slow rotators. This 
effect appears to be 1.3 to 3 times magnified in the DR-uncorrected VI and BI CMDs, ac- 
cordingly, which can be attributed to a likely anticorrelation between the rotational blueshift 
and DR. 

To validate these findings, we explored the likelihood of a chance blueshift in our sample. 
We constructed a synthetic cluster with all physical parameters close to those of our selection 
of Tr 20 members, particularly, the total number of stars on the main sequence and the 
pattern of DR. The only variable was a level of reddening "noise" added to the grid of median 
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DR (Figure 3). From a 1000 realizations of a synthetic cluster, each time we randomly chose 
168 stars and calculated the median V — I color difference between the "fast" and "slow" 
rotators. In the absence of stellar rotation, the standard deviation of these median differences 
is e{V — I) = 0.005 and it is only weakly dependent on the level of the applied DR-noise. 
Thus, in a statistical sense the estimated color offset due to the stellar rotation is a ~ 2(j 



result and agrees qualitatively with the theoretical predictions by lGirardi et al.l (120111 ). 



Near the MSTO, a 6{V — I) = —0.01 color offset for the 1.3 Gyr isochrone corresponds 
to an age shift of ~20 Myr. Considering the median difference of ~60 km s~^ between the 
slow and fast rotators and extrapolating it to the entire range of vsini (180 km s~^), we 
obtain a tentative "dispersion" in isochron e ages of ^60 Myr. Th is is only ~15% of the 



deduced average spread of ages reported by iGoudfrooij et al.l ( 120111 ). We note that contam- 



ination of the slow-rotator sample by field stars and fast rotators at high inclination angles 
masquerading as slow rotators, can reduce the true color dispersion by ~0.01 mag. Even 
accounting for this, the stellar rotation appears to have a marginal effect on the morphology 
of MSTO. 



5. Conclusions 

This study of the Galactic open cluster Tr 20 provides strong evidence that DR can 
mimic the appearance of multiple populations in star clusters. Despite some earlier indica- 
tions, Tr 20 appears to contain a single coeval population of stars. Positions of the upper 
main-sequence stars in the CMD of Tr 20 are marginally dependent on their projected ro- 
tational velocities and in a specific way. Our data show that high vsini stars are slightly 
blueshifted. Stellar rotation appears to play merely a minor role in the observed broadening 
of MSTO of this intermediate-age open cluster. 

We thank Nate Bastian and the anonymous referee for meritorious suggestions. This 
work has been supported in part by the NSF Grant AST 09-08114 to JHU (LP.). 
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Fig. 1. — Upper panel: radial velocities as a function of f sini. The large dots indicate our 
selection of probable kinematic cluster members of Tr 20. The widening range of cluster 
Vj- is due to the growing uncertainties at higher vsini. Lower panel: a histogram of the 
main-sequence probable cluster members. A sample used to examine the effects of stellar 
rotation (see Section 4) is marked by shaded area. 
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Fig. 2. — Color-magnitude diagram of Tr 20. Left panel: the small dots show all stars within 
r < 8' from the cluster center. The circles indicate our selection of targets for VLT/FLAMES 
spectroscopy. Middle panel: the larger dots are \4-selected probable cluster members. The 
arrow shows the direction and a half of the adopted amount of reddening, E(y — I) = 0.50. 
The location of the upper main sequence at this lowest reddening is indicated by the curve. 
The presumable field stars are deleted at \^ < 17. Right panel: all stars shown in the middle 
panel but corrected for differential reddening. 
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Fig. 3. — Spatial distribution of differential reddening (DR) in Tr 20. Upper panel: estimated 
DR using the main-sequence stars with vsini < 95 km s~^. The size of filled circles is 
proportional to the amount of DR, ranging from 0.0 to ~0.13 mag in E{V — I). A window 
of lower reddening is clearly visible at X ~ —1', Y ~ +3'. Lower panel: a grid of smoothed 
DR. In the corners of this grid the estimates of DR are extrapolated and uncertain. 
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Fig. 4. — Color-magnitude diagram of Tr 20 and a fit with theoretical isochrones. Left panel: 
same as Figure 2 (middle panel). Provisional fit with a Padova 1.3 Gyr isochrone targets the 
cluster's "blue" main sequence. A poor fit to the features of RGB is evident. Middle panel: 
same as Figure 2 (right panel) and the overplotted Padova 1.3 Gyr isochrone. The tightness 
of a fit has improved substantially. There is evidence that the morphology of the red clump 
is complex. Right panel: a fit with a Dartmouth 1.5 Gyr isochrone. The reasoning of an 
older age is given in the text. In all three cases the fit parameters (see Section 3.2) are kept 
unchanged. 
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Fig. 5. — Color distribution of the upper main-sequence stars. The sample of 168 stars 
with V < 16.3 is divided at t>sini = 90.5 km s~^ in two equal halves of slow (dark curve) 
and fast (red curve) rotators. The dotted curves represent the normalized and smoothed 
histograms after correcting the abscissa for differential reddening, while continuous curves are 
constructed with uncorrected data. In both cases, the fast rotators are slightly blueshifted. 



